Phase Diagram of the Proton-Neutron Interacting Boson Model 
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We study the phase diagram of the proton-neutron interacting boson model (IBM-2) with special 
emphasis on the phase transitions leading to triaxial phases. The existence of a new critical point 
between spherical and triaxial shapes is reported. 
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Quantum phase transitions (QPT) have become a 
subject of great interest in the study of several quan- 
tum many-body systems in condensed matter, quantum 
optics, ultra-cold quantum gases, and nuclear physics. 
QPT are structural changes taking place at zero temper- 
ature as a function of a control parameter (for a recent 
review see 0). Examples of control parameters are the 
magnetic field in spin systems, quantum Hall systems, 
and ultra-cold gases close to a Feshbach resonance, or 
the hole-doping in cuprate superconductors. 

The atomic nucleus is a finite system composed of N 
neutrons and Z protons (Z+NwlOO). Though strictly 
speaking QPT take place for large systems in the ther- 
modynamic limit, finite nuclei can show the precursors 
of a phase transition for some particular values of N and 
Z. In these cases, one finds specific patterns in the low 
energy spectrum revealing the strong quantum fluctu- 
ations responsible for the phase transition Q. Recently 
the concept of critical point symmetry has been proposed 
by Iachello and applied to atomic nuclei. First, the tran- 
sition from spherical to deformed 7-unstable shapes was 
studied and the corresponding critical point called E(5) 
0- Since then, the interest in nuclear shape-phase tran- 
sitions has been constantly growing. The characteristics 
of the critical point in the phase transition from spherical 
to axially deformed nuclei, called X(5), were presented in 
Ref. Q. More recently, the critical point in the phase 
transition from axially deformed to triaxial nuclei, called 
Y(5), has been analysed [fj. In all these cases, critical 
points are defined in the context of the collective Bohr 
Hamiltonian Q . Using some simplifying approximations 
precise parameter-free predictions for several observables 
are obtained. This allows to identify nuclei at the critical 
points looking at spectroscopic properties. Indeed, some 
experimental candidates to critical nuclei have already 
been proposed 0,0- 

The collective Bohr Hamiltonian, underlying this ap- 
proach to critical point symmetries, is closely related to 
the Interacting Boson Model (IBM) 9]. The simplest 
version of the IBM is called IBM-1 since in it no ex- 
plicit distinction is made between protons and neutrons. 
In IBM-1 there are three dynamical symmetries: SU(5), 
0(6), and SU(3). These correspond to well defined nu- 
clear shapes: spherical, deformed 7-unstable, and prolate 



axial deformed, respectively. The structure of the IBM-1 
Hamiltonian allows to study systematically the transition 
from one shape to another. There were some pioneering 
works along these lines in the 80's 0, 0, 0, but it 
has been the recent introduction of the concept of criti- 
cal point symmetry that has recalled the attention of the 
community to the topic of quantum phase transitions in 
nuclei. The phase diagram of the IBM-1 has been stud- 
ied from several points of view [3 13 13 13 13 El 
The three different phases are separated by lines of first 
order phase transition, with a singular point in the tran- 
sition from spherical to deformed 7-unstable shape that 
is second order. In the usual IBM-1 no triaxial shapes 
appear. These can only be stabilised with the inclusion 
of specific three body forces. A more natural way to 
generate triaxial deformations is by explicitly taking into 
account the proton-neutron degree of freedom with the 
more realistic IBM-2 [3]. 

In this letter we will study the phase diagram of the 
IBM-2 using a simplified Hamiltonian that keeps all the 
main ingredients of the most general one. This is the 
Consistent-Q IBM-2 Hamiltonian [HI 
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- Xn d\,d K and N is the to- 
tal number of bosons, which is equal to the number of 
valence proton plus neutron pairs. The IBM phase di- 
agram studied up to now corresponds to the selection 
Xtt = Xv which produces cither spherical, axial or 7 in- 
dependent shapes. We will extend the previous works on 
IBM phase transitions by exploring the transitions from 
axial to triaxial shapes within the mean field or intrin- 
sic state formalism. The trial wave function is the most 
general proton-neutron boson condensate 0, IT^ . |20| . 
\g) = \N^,N v ,f3 v , 7^,^,7^,0) 
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where n = Tr,i> and i?3(f2) is the three dimensional ro- 
tation operator with f2 fixing the relative orientation 
(Euler angles) between the proton and neutron conden- 
sates. Nn and N v are the numbers of valence proton and 
neutron pairs, respectively. The equilibrium values of 
the structure parameters (fln,^, j3 v ^ v , Q) and the en- 
ergy of the system for given values of the control pa- 
rameters in the Hamiltonian (x, XiriXv) can be obtained 
by minimising the expectation value of the Hamiltonian 
||TJ in the intrinsic state J2J : 6{g\H\g) = 0. Although 

x=0 X=-\=^F>/2 
SU*( 3 ) 
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FIG. 1: Pictorial representation of the IBM-2 parameter space 
with a dynamical symmetry in each of the four vertices. 

there is an explicit dependence of the energy on the Eu- 
ler angles, it has been shown [2(j that oblique config- 
urations (relative orientation angles different from the 
aligned Q = or the perpendicular Q — tt/2 ) require 
a repulsive hexadecapole n — v interaction. Therefore, 
since our Hamiltonian has no hexadecapole terms, 
we do not expect oblique configurations. We can then 
safely assume that any arbitrary local minimum will 
have SI = 0, 7^ = 7„ = 0° (or equal to 60°) for the 
aligned configurations or £1 = 0, 7^ = 0°, 7„ = 60° (or 
7^ = 60°, 7^ — for the perpendicular configurations. 
In both cases O = and the rotation operator disappears 
from the intrinsic state J2J- In that situation, the energy 
per boson in the limit N„ , N v — > 00 reduces to 
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where we have used the notation Qq(k) = Qo(A 7,x)k 
and Q 2 (k) = Q- 2 {n) = Q2(P,H,x)k = Q-2(J3,7,x)k 
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FIG. 2: Transition from 5(7(3) to SU* (3): x = 0, x> = 
—\/7/2 and Xv varies from — V7/2 to +y/7/2. In the panels 
are plotted the energy of the ground state in arbitrary units, 
and the variation of the shape parameters (5 V (dimensionless) 
and 7„ (degrees). 



As a natural extension of the Casten triangle for IBM- 
1 01 > the geometrical representation of the IBM-2 is 
a pyramid with the new triaxial dynamical symmetry 
SU*(3) j2jj in the upper vertex. Fig. ^ shows a pictorial 
representation of the IBM-2 parameter space. Any point 
in this space is obtained with the following transforma- 
tion to polar coordinates (see Fig. HJ. 
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We have explored the IBM-2 parameter space of Hamil- 
tonian and present here a selected set of calculations 
in order to establish the IBM-2 phase diagram (a more 
detailed presentation will be given in a forthcoming pub- 
lication). We have not found traces of phase transitions 
in the transition from 0(6) to SU*(3) in a parallel way 
as the already known transition from 0(6) to SU(3) in 
IBM 1. The 0(6) symmetry is in fact very unstable 
against small perturbations driving the system out of the 
dynamical symmetry either to axial deformed or to triax- 
ial shapes depending on the interaction. The 0(6) sym- 
metry itself has been proposed as a critical dynamical 
symmetry [22j. 

In Fig. we show the transition SU(3) -> SU*(3) 
through the edge plotted in Fig. 1. Along this line 
x = and x-n = — V"/2 are fixed. The relevant con- 
trol parameter is Xv varying from — v7/2 (equal and 
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aligned quadrupole prolate shapes for protons and neu- 
trons) to y/7 /2 (quadrupole prolate shape for protons and 
quadrupole oblate shape for neutrons with perpendicular 
axis of symmetry |2l|). In Fig. [21 we present the results 
for the ground state energy (in arbitrary units) and the 
shape parameters (/3„,7„). The resulting proton param- 
eters are (3^ = v2 and 7^ = for all values of the control 
parameter In the limit \ v — —y7/2 we recover the 
results known from IBM-1: /3 U — \/2 and j u = 0. In 
the opposite limit Xv — V< /2 the results known from 
Ref. [2l| are obtained: /3„ — V2 and 7^ = 60°. Around 
Xv = 0.4035 a clear shape phase transition is observed 
changing the system from axial {xu < 0.4035) to triaxial 
(Xu > 0.4035). We will call this point "y". Note that 
in this phase transition the order parameter is 7„ that 
changes from in the symmetric phase to a finite value 
in the non-symmetric phase |23|. We have minimised 
the energy following two inverse paths looking for pos- 
sible coexistence of minima and the corresponding spin- 
odal and antispinodal points. Both calculations give ex- 
actly the same results. This means that spinodal, critical 
and antispinodal points all converge to a single point [2J . 
Therefore, the transition from SU(3) to SU*(3) is second 
order. 
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FIG. 3: Same as Fig. but for the transition from U(5) to 
SU*(3): \t ~ ~Xv = — vi /2 and x varies from (triaxial) 
to 1 (spherical). 



Fig-EHshows the transition from U(5) — > SU*(3) through 
the corresponding edge in Fig. ^ Along this edge x-n = 
~Xv = — v7 /2 are fixed and the relevant control param- 
eter is x changing from 1 (spherical) to (triaxial). The 
values of /3 ff and (3 U are always equal at the energy min- 
imum, while 7 T and j u are symmetric with respect to 
7 = 30° axis. In the different panels, the energy, and 
the values of (3 and 7 for proton and neutron shapes are 
presented. For x = 1, fl^ = /3„ = implying a spheri- 
cal shape. For x = we recover the SU*{3) case with 
Pit = flu = v2j and 7^ = and 7„ = 60°. A phase tran- 



sition at x = 0.8 is observed. We will call this point "x*" . 
As in the preceding case, we have performed two sets of 
calculations following inverse paths to determine the or- 
der of the transition and again we have found no region 
of coexistence, converging at the same place, spinodal, 
critical, and antispinodal points. Note that in this case 
the order parameter is /3 n = {3 V , as well as 7^ = 60 — 7„ . 
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FIG. 4: Same as Fig. |2]but for a generic transition from U(5) 
to a triaxial shape. The structure parameters are x-n = —1-2, 
Xv = 0.5 and the control parameter x varies from to 1. 



In Fig. 01 we present the study of a generic transition from 
U(5) (spherical) to a triaxial shape through a trajectory 
within the IBM-2 pyramid. In particular, we have se- 
lected the trajectory defined by Xn = —1-2 and Xv = 0.5, 
using x as the control parameter varying from 1 to 0. 
The ground state energy, and the values for the /?'s and 
the 7's are plotted. Two phase transitions are observed 
at different values of x. Starting from x = 1 (spherical 
system), a first transition to axial deformed shape is ob- 
served at x ~ 0.8. At this point, the values of (3„ and 
/3 V depart from zero but 7^ and 7„ are zero indicating 
a deformed axial symmetry. [3^ and /?„ play the role of 
order parameters in this phase transition. For a value 
of x w 0.48 a second phase transition is observed. The 
values of f3 w and /3 V are different from zero in both sides 
changing smoothly along the transition. The angular pa- 
rameters 7 jump from zero to finite values indicating a 
transition from an axial shape to a triaxial shape. There- 
fore 7,r and 7„ are the order parameters. The different 
values for the shape parameters for protons and neutrons 
are due to the selection of the structure parameters x™ 
and Xv for this trajectory. As in preceding cases we have 
performed two sets of calculations following inverse paths 
to determine the order of the phase transition. The tran- 
sition at x w 0.8 can be analysed looking at the behaviour 
of the ground state energy (inset). Where the full line 
corresponds to a forward calculation, starting at x = 0, 
increasing x, and the dashed line to a backward calcu- 
lation, starting at x = 1, decreasing x. The inset shows 
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FIG. 5: Schematic phase diagram for IBM-2. S stands for 
spherical, A for axial and T for triaxial phases. The critical 
points "x*" and "y" studied here and those already known for 
the IBM-1 phase diagram, "x", "e", and "0(6)", are marked 
with dots. 

that there are two minima competing, one spherical and 
one deformed. If the system comes from the spherical 
region it keeps spherical for a while even there is another 
deformed minimum with slightly lower energy. On the 
other side, if the system comes from the deformed re- 
gion it keeps deformed (look at the small peak in the full 
line in the inset at x=0.802) although another spherical 
minimum have slightly lower energy. This coexistence of 
deformed and spherical minima in a small region around 
x = 0.8 is the signature for a first order phase transi- 
tion. The phase transition at x m 0.48 has been studied 
with forward and backward calculations. We have not 
found any coexistence region. The antispinodal, critical 
and spinodal points come together to a single point as 



corresponds to a second order phase transition. 



We have explored the parameter space of the IBM-2 
Hamiltonian in Fig. 2] The resulting phase diagram 
of the proton-neutron IBM as described by the Hamil- 
tonian m is depicted in Fig. [S] There are three well 
defined phases: spherical, axially deformed (prolate in 
the schematic presentation of Fig. [SJ) and triaxial. The 
critical surface separating spherical and axially deformed 
shapes (e-x*-x-e) is first order, while the surface separat- 
ing axially deformed and triaxial shapes (e-0(6)-y-x*-e) 
is second order, including the common line between both 
surfaces (e-x*). We have checked that in all the cases 
discussed in which the transition is second order the be- 
haviour of the corresponding order parameter near the 
critical point is consistent with a critical exponent 1 /2 as 
given by the Landau theory [23] . We would like to stress 
that the critical surface separating spherical and axially 
deformed nuclei is almost a sphere with a radius equal to 
p = 0.2 and centered in U(5). The straight line plotted 
inside the figure gives an idea of the trajectory followed 
by the transition discussed in Fig. 0] We would like to 
emphasise that we have found a new critical point (x*) at 
the phase transition changing directly from spherical to 
triaxial shapes. We are currently studying the spectro- 
scopic properties of this critical point. The results will 
be presented elsewhere. Finally, this scheme of analysis 
can be easily extended to positive values of \-n to obtain 
the dynamical symmetry limits SU{?>) and ST/* (3). 
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